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For Kraichnan's problem of passive scalar advection by a velocity field delta-correlated in time, 
the limit of large space dimensionality d 3> 1 is considered. Scaling exponents of the scalar field are 
analytically found to be (2n = nC,2 — 2(2 — Qi)n{n — l)/d, while those of the dissipation field are 
fi„ = —2(2 — (,2)n(n — l)/d for orders n <C d. The refined similarity hypothesis C^2 n — n (,2 + is 
thus established by a straightforward calculation for the case considered. 
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It is likely that Kraichnan's model of white-advected 
passive scalar jjj will become a paradigm in theoretical 
studies of intermittency and anomalous scaling in turbu- 
lence . This is because any simultaneous correlation 
function of a scalar satisfies a closed linear differential 
equation so that all common hypotheses about intermit- 
tency could, in principle, be verified by direct calcula- 
tion. In an isotropic turbulence, the 2n-point correlation 
function depends on n(2n — 1) distances, which makes 
direct solution of the respective partial differential equa- 
tion quite difficult in a general case. The fourth-order 
correlation function has been calculated recently in two 
limiting cases: i) large space dimensionality d 3> 1 |0] and 
ii) almost smooth scalar field 2 — £2 <C 1 Q . By £ 2n we 
designate the leading scaling exponent of the structure 
function: S-2 n {r 12) = — 2 ) 2n ) oc r£ 2 n + subleading 
terms. If (2™ 7^ nCz then the scaling is called anoma- 
lous. Anomalous scaling is the way intermittency man- 
ifests itself in developed turbulence: the degree of non- 
Gaussianity, that may be characterized by the dimen- 
sionless ratios S^n/S^, depend on scale. We shall see 
below that the anomalous dimensions A 2 „ = n£ 2 — Qin 
are positive which means that the smaller the scale of 
fluctuations the more non-Gaussian the statistics is. 

In this Letter, we use the formalism developed in Ref. 
3 to calculate high-order correlation functions assuming 
l/d to be the smallest parameter in the problem. Follow- 
ing , we demonstrate how the anomalous part of the 
solution appears as a zero mode with the form indepen- 
dent of the pumping. Those zero modes exploit the in- 
terchange symmetry between the distances, the number 
of zero modes and anomalous exponents thus increases 
with the order of the correlation function. 

The advection of a passive scalar field 8(t, r) by an in- 
compressible turbulent flow is governed by the equation 



kA)6 = 4>, 



(1) 



V a u a — being assumed. The external velocity u(i, r) 
and the source <f>(t, r) are independent random functions 
of t and r, both Gaussian and (^-correlated in time. 
Their spatial characteristics are different. The source 
is spatially correlated on a scale L i.e. the pair corre- 
lation function ri)0(t 2 , r 2 )) = 8{ti - t 2 )x{ r i2) as 
a function of the argument r\ 2 = |fi — r 2 | decays on 



the scale L. The value x(0) = P is the production 
rate of 6 2 . The velocity field is multi-scale in space 
with a power spectrum. The pair correlation function 
(«"(ii,ri)/(i 2 ,r 2 )) = S(h - t 2 )[V 6 a P - /C^(r 12 )] is 
expressed via the so-called eddy diffusivity 
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where < 7 < 2 and isotropy is assumed. 

Considering steady state and averaging (|l|) over the 
statistics of u and 4> one gets the closed balance 

equation for the simultaneous correlation function of the 
scalar Fi... 2n = F(r x , . . . ,r 2 „) = (0(ri), . . . ,0(r 2n )): 

-£Fi...2n=Fi---2n-2X2n-i,2n +permutations. (2) 

The operator t = £V : /C Q/3 (r y -)Vf V^/2 + nJ2 A,: de- 
scribes both turbulent and molecular diffusion, it may 
be rewritten in terms of relative distances ra fl3|] : 
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Here, the summation is performed over n(2n— 1) indepen- 
dent distances (for d > 2n — 2) with subscripts satisfying 
the conditions i^j and m ^ i,j , n ^ the diffusion 
scale r 2 , 1 = 2K(2—j)/(D(d — 1) has been introduced. 
We consider the convective interval L 3> 3> where 
the operator C is scale invariant. Taking n = 1, it is easy 
to find the pair correlation functions of the scalar |]|| : 



F 12 (r) = P 
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We see that (2 — 7- The scaling exponent of C is —7; 
the solution of (||) may thus be presented in the form 
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F = Ff orc + Z, where we separated the so-called "forced" 
part of the solution (with the scaling exponent C271-1 + 7 
prescribed by the rhs) from the zero mode Z that may 
have a different scaling. It has been recognized inde- 
pendently by the authors of ]3|-|| that they are the zero 
modes of the operator £ that are responsible for the 
anomalous scaling. We shall demonstrate below that the 
factor (L/r) A2n appears in the zero mode while the fac- 
tor (L/r) A2n - 2 appears in the forced term, A2„ being the 
anomalous dimension. The zero mode turns out to be 
dominant. It has been demonstrated in H that does 
not appear in the leading terms of F as long as at least 
some of the distances are in the convective interval. 
Assuming that to be the case, we omit the diffusive parts 
of £. We shall account for them later while considering 
the correlation functions of the scalar derivatives. 

Let us consider now the case of large space dimen- 
sionality where the anomalous dimensions can be cal- 
culated analytically. It is seen from (^) that the level 
of scalar fluctuations necessary to provide for a given 
flux P decreases as d increases. It may be shown that 
F\...2n oc P n /[d(d— 1)]™. Considering large d and assum- 
ing that the flux is determined by the pumping (and it 
is thus d-independent) we shall develop the perturbation 
theory for the quantities dr n Fi,..2 n which have finite lim- 
its at d — > 00. We shall show below that despite the small 
level of fluctuations at large d, the statistics of the scalar 
is substantially non-Gaussian at small scales. However, 
since the anomalous dimensions are small, there exists a 
wide interval of scales where the correlation functions are 
close to their Gaussian values and can be calculated by 
perturbation theory to obtain A.2n at the leading order 
in 1/d. As d — > 00, the main part of £ is the operator of 
the first order £ =(d 2 D/ (2 — 7)) Yfij^ 7 ^r i:j ■ Since £ is of 
the second order, one may wonder if the 1/d perturbation 
theory is regular. It has been demonstrated in j!| that it 
is so by considering the perturbation theory that starts 
from the bare operator £' cx r ij~ d dr ij 'r]^ d d rij ; this 
alternative theory gives the same answer for the anoma- 



lous exponent. The difference between the perturbation 
theories with £q and £' Q is in the representation for the 
forced term (the second approach corresponds to some re- 
summation). Since we are interested here in the anoma- 
lous scaling which is given by the zero modes that are 
determined by their behavior at r = 0, there is no quali- 
tative difference between £q and £' and the perturbation 
theory for the anomalous exponent is regular. 

The zeroth term in the perturbation series for F\...2 n is 
given by a Gaussian reducible expression. Our aim is to 
iterate it once by applying the operator £q 1 (£— £q). The 
parameter of the expansion is n/^/d, assumed to be small. 
In the first correction to be thus found, the logarithmic 
terms ln(L/ry ) are of interest because they appear at 
expanding the anomalous scaling factors (L/rij) A2n over 
— see [^) for the details. As we shall see, there are 
different zero modes with different anomalous dimensions 
^2n.i for any given n. Of course, only the largest &.2n,i 
contributes in the limit L/r — > 00. However, in the region 
^2n,i lu(^/ r ij) 1, where we carry out our calculations, 
all A's contribute logarithmic terms. We thus have a de- 
generate perturbation theory (all zero modes Z2n,i have 
the same exponent in the zeroth approximation) and 
should therefore apply the operator £§ x {£ — £q) on a 
vector of zero modes of £0 and then single out the terms 
having logarithms. It is clear that logarithms may appear 
only multiplied by a zero mode of £q. We thus obtain 
the matrix of the operator £q 1 (£ — £q) m the represen- 
tation of Z2n,i- The eigenvalues of that matrix are the 
anomalous exponents A2 n ,i at the leading order in 

Let us describe how the matrix is generated. The 
most convenient classification of the zero modes is as fol- 
lows: Z2n.i is the polynomial in x — r 1 of order n which 
may be separated into a symmetrical sum of polynomials, 
each involving distances between i points. For example, 
there are two zero modes for the fourth-order correlator : 

= J2( x ij ~ x ki) 2 and Z 4j3 = J2( x ij ~ x jk) 2 f|- The 
first-order logarithmic correction is calculated by the rule 



-£ {£-£ )x t3 x kl - 2d m[rjj 2(x ik -x a -x jk +x jl )*, ./j^/i, (5) 



which gives the matrix || 

Here, A 4 , 4 = 4(2 ~i)/d and A 4 , 3 = -(4- 7 2 )/2d arc the 
eigenvalues. Before describing the structure of the matrix 
at higher n, let us note that the order n and the number 
of points i are not enough to specify the zero mode for 
n > 3 and 4 < i < 2n — 2 due to the possibility of differ- 
ent topological configurations (we enumerate them by j). 
For example, at n = 3, i = 4, there are two different zero 



modes: one involves distances rki,ri m ,r mn and another 
tm ■, fin , fim . The total number of zero modes grows with 
n faster than factorially due both to the growth of the 
number of possible functional forms and the appearance 
of new topologically different configurations. What is im- 
portant to know is that the operator acting on the zero 
mode with a given i produces only modes with i' < i. 
The analysis of eigenvalues is thus reduced to the consid- 
eration of the blocks with a given i. The first mode Zm^n 
contains the monomial 0:1,2 ■ ■ ■ %n-i,n that cannot be ob- 
tained by da) from other modes, so that the first element 
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is n(n — l)A4 ) 4/2 and the remaining elements of the first 
column are zero. Considering -H2n,2n-ii one realizes that 
the second diagonal element is (n— l){n— 2)A 4i 4/2 + A4 j3 
and all the lower elements in the second column are 
zero. Then the 3x3 block follows which corresponds 

to Z2n,2n-2J- 



. /(n-l)(n-4)+4< ? 
A!l ' 2(3 -n) r 




-2 
-3n-l-6g l + Aq 
9 + l2q (ra-2)(ra-5), 



(7) 



with q = A4 j 3/A4.4, all the elements below the block are 
zero. The next block is for i = 2n — 3, it is 7 x 7 for 
n > 6, the sizes of the blocks grow as one approaches 
the center of the matrix, then they decrease and even- 
tually we come to the 3x3 block due to Z2n,2,j and 
a single value at the lower right corner. By virtue of 
(||), all entries of the matrix at arbitrary n may be ex- 
pressed via those A4 4 and A4 3 and combinatorial fac- 
tors. Note that the classification of all eigenvalues and 
eigenvectors is very important because it carries informa- 
tion about the algebraic structure and underlying sym- 
metry that governs our set of correlation functions. We 
postpone the general classification until further detailed 
publications. Fortunately enough, the mode that gives 
the largest anomalous dimension and the structure func- 
tions of the dissipation field is separated so that it can be 
found without finding the whole set of A2 n ,i- We notice 
that since A4 j4 > A4 ; 3 and the largest combinatorial fac- 
tor in front of A4 t 4 appears in the first element then it is 
plausible to assume that Z 2n2n gives the largest eigen- 
value A2„.2n = n(n — l)A^i/2. One can directly check 
that A2 ni 2n-i and all eigenvalues A2 n ,2n-2 of (Q) are less 
than A2 n> 2ri for any n. For an arbitrary block, the valid- 
ity of the assumption may be established asymptotically 
for 71 ^ 1 (yet, of course, n <C jd) when all eigenval- 
ues are 71^4,4/2 + 0(n). For n = 2,3,4, we found all 
the eigenvalues using Mathematica; the largest is always 
A 2 „. 2n . We thus conclude that for A2 n) 2n ln(L/r) 1, 
the main contribution to the correlation function is given 
by the respective zero mode with the scaling exponent 
&n = nC,2 — 2(2 — C,2)n(n — l)/d. In particular, 



((01 - e 2 f n ) ~ r »7(£/ r )2«(«-l)(2-7)Ai_ 



(8) 



It agrees with wm where £4 has been calculated. Note 
that both (|) and the results of [||J^] differ from £ 2n sug- 
gested in p I . In our opinion, that means that the closure 
implemented in 0] is not exact at the limits considered 
(n <C jd and 2 - 7 <C 1). 

To find the correlation functions of the scalar deriva- 
tives, one should consider some distances ry as going to 
zero. While some distance passes the diffusion scale the 
dependence on that distance changes. To describe that, 
we include the diffusion operator into C As a result, the 
diffusion scale r d appears in the correlation functions. 
The form of r d dependence could be readily established 



for an arbitrary n, d, 7 by using a straightforward per- 
turbation expansion in the ratio between small and large 
distances (see Sect. Ill of ||). The overall scaling of 
F\...2n at all the distances inside the convective interval 
is assumed now to be known (A2„ = A2„ i 2n for d 1) 



F\...2n ~ C 2n R 



717— A 2 



-A,, 



R-»r d . (9) 



Now, let us consider one distance, say ri2, to be much 
smaller than the others: p = r\ 2 <C rij ~ R. At zero 
order in p, F 1X3 ... 2n « G(R) ~ #™7-A2„£A 2 „ || Thc 
leading isotropic correction satisfies the equation 



C p SF 2n (R, p) = C R G(R) - $(R) 



(10) 



where $(i?) - #("-i)7-A 2 „- 2 £A 2n _ 2 ig the ma j or 
dependent) term of the rhs of (JsJ) , is the major term 
of the operator C, and C p is the perturbation operator 
t p = /C a/3 (p)V£V^ + 2kA p , . The solution has the form 



SF 2n (R,p) _ i? («-i)7-A 2 „ L A 2 „ f 

Jo 



rdr 



2r 2 " 7 



+ r 2 -7 



(11) 



In deriving (|1 1[) , it has been implied that A 2 „ > A2 n _2- 
At p <C r d , the isotropic correction is analytic in p: 



,7-2 p 2 



Now we are 



5F 2n {R lP ) ~ ij(n-i) 7 -A 2 „ L A 2 
ready to differentiate it with respect to p. in particular, 
to calculate the correlation functions that involve the dis- 
sipation field e = k[V9] 2 : 



1P3 



•0 2n » ^R (n -^{L/R) 



A 2 



(12) 



Since e oc K oc r 2 , 1 then r d disappears. The diffusion 
scale appears only at the anisotropic terms proportional 
to angular harmonics with respect to the angle between 
R and p. This is because of the zero modes of the oper- 
ator C p associated with the angular harmonics that are 



expressed via the Jacobi polynomials P, 



2 k 



z 2fe ( P ,R) = p 2 ( r ) 



"np" 






. P _ 





p 5 ", 

„2fc _«k-2fc 



p > r d , 
p < r d ; 



5 k=2\ l-d+\j{d- 1 f + 



8/c(d+l-7)(2fc + d-2) 



where v = (d — 3)/2, n = R/R. The 2fc-th angular cor- 
rection to fi 1 3,... 2n is thus 

S^F 2n (p^) ~ Z 2k {p,R){L/R) A2n R n ^ s \ (13) 

where the leading scaling behavior of the p-independent 
prefactor is restored to give a correct dimensionality. The 
overall (L-dependent) scaling is assumed to be known. 
If 717 > 5k and p is in the diffusive interval, then the 
estimate (|l3|) is true only at R being small enough: 
(R/L) A ^(r d /R) n ^~ s " < 1 — see @ for the details. 
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Differentiating 8^ 2k 'F with respect to p one can es- 
tablish the scaling of the correlation functions that in- 
volve high-order traceless tensors of derivatives £} k ^ = 
K (2fe-7)/(2-7)^ =oai [( n v) fc -^V ai •••V a <0] 2 ! m is a re- 
spective coefficient of f 2 fc (%) expansion in the series in 
x. One can also generalize (II4 ) in another way, consid- 
ering extra fusion of another one or more pairs of points. 
The results show that the fusion procedures for an arbi- 
trary number of different pairs of points commute with 
each other. In other words, using the language of the 
so-called operator algebra P~|fl"f the validity of which for 
turbulence was argued in [ |12[|13[| , we can introduce r d - 
related dimensionality of the fields: the scalar field 9 and 
dissipation field e have dimensionality 0, while £}- k ^ has 
dimensionality 8k — 7- Thus, to find the ultraviolet (r^- 
related) dimensionality of a composite multiplicative field 
one should sum the dimensionalities of the multipliers: 



(k)]c k 



oc r d 



The above fusion rules describe the scaling of the fluc- 
tuations of the enumerated primary fields only. Partic- 
ularly, the rules give the scaling of correlation functions 
with e at different points but not with the higher powers 
of the dissipation field, e™. The difference stems from the 
fact that to calculate the correlation functions contain- 
ing e n we should fuse a group of 2n (rather than two) 
points in an initial 9 correlator. To calculate a correction 
that produces a nonzero contribution into the desirable 
correlator of e n , one should compare the leading forced 
correction with a zero mode of the 2n-point operator, but 
not with that of two-point one C p . Such zero modes do 
have anomalous scaling, as we have learned above, they 
produce the major contribution to correlation functions 
of e". The ultraviolet anomalous scaling of such objects 
is thus related to the infrared anomalous scaling of the 
passive scalar itself. For ex amp le, by a direct application 
of the above procedure ( |n)||L2| ), we get 

N A 2 „+A 2 



l (L/r d ) 



(14) 



The dissipation field is thus highly intermittent, the 
single-point means (e n ) ~ {e) n (L /rd) A2n grow unlimited 
when diffusivity decreases. Statistics at the convective 
interval is better related to local average e r over the ball 
with the radius r. Since spatial integration and time av- 
erage commute then our knowledge of (ei ■ ■ ■ e n ) with all 
distances in the convective interval allows one to obtain 
by spatial integration (which converges if A4 < d) : 

((e r ) n ) ~ (e>"(r/L)"», /x„ = C 2 „ - n( 2 = -A 2n . (15) 

That relation presents a version of the refined similarity 
hypothesis ]l5|-p^| valid in our case. Indeed, scalar field 
and dissipation field are related: e r ~ (89 r ) 2 /r r where 
r r ~ Dr 1 is a transfer time. For n — 2, ( JTa ) has been 
established in |||. 



The first (n^jd) moments of the locally averaged dis- 
sipation have fx n — —n(n— l)A 4j4 /2 so they are described 
by /03-normal statistics [n4j : 



(ln[e r /(e)] -A 4 , 4 ln[i/r]/2) i 



V{e r ) ~ exp - 



2A M ln[L/r 



To conclude, the law (g) qualitatively corresponds to 
the observed behavior of £ n [p^|-[20| . Quantitatively, we 
cannot use (||) for d = 3, 7 = 2/3 (which would give an 
overestimation of the anomalous exponents) because the 
validity condition of our theory n C 7<i will be violated 
already at n — 2. Note that quadratic dependence of ( n 
and [i n on n (and lognormality) is violated when n ~ jd, 
while the similarity relation ( |l5| ) is true for any n, 7, d. 

This work is a part of the extensive program on study- 
ing anomalous scaling in turbulence undertaken together 
with E. Balkovsky, I. Kolokolov and V. Lebedev. We 
are grateful to them for numerous discussions. We are 
indebted to G. Eyink, U. Frisch, R. Kraichnan and the 
anonymous referee for useful remarks. The work was sup- 
ported by the Clore Foundation (M.C.) and by the Rashi 
Foundation (G.F.). 
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